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APPROXIMATION  OF  THE  OPTIMAL  COMPENSATOR  FOR 
A  LARGE  SPACE  STRUCTURE1*2 


Michael  K.  Mackay 

University  of  California,  Los  Angeles 


ABSTRACT 


'This  paper  considers  the  approximation  of  the  optimal  compensator  for  a 
Large  Space  Structure.  The  compensator  is  based  upon  a  solution  to  the  Linear 
Stochastic  Quadratic  Regulator  problem.  Colocation  of  sensors  and  actuators  is 
assumed.  A  small  gain  analytical  solution  for  the  optimal  compensator  is  obtained 
for  a  single  input/single  output  system,  i.e*,  certain  terms  in  the  compensator 
can  be  neglected  for  sufficiently  small  gain.  The  compensator  is  calculated  in 
terms  of  the  kernal  to  a  Vol terra  integral  operator  using  a  Neumann  series.  The 
calculation  of  the  compensator  is  based"  upon  the  b^semlgroup  for  the  Infinite 
dimensional  system.  A  finite  dimensional  approximation  of  the  compensator  is, 
therefore,  obtained  through  analysis  of  the  infinite  dimensional  compensator 
which  is  a  compact  operator. 


'•A 


1.0  INTRODUCTION 

One  of  the  distinguishing  properties  of  a  Large  Space  Structure  (LSS)  is 
that  it  is  a  distributed  parameter  system  and,  hence,  an  infinite  dimensional 
mathematical  model  is  required  for  its  description.  In  most  applications, 
active  control  of  shape,  attitude,  and  structural  vibrations  will  be  necessary. 

For  such  problems,  formulation  of  the  control  problem  as  a  steady-state  (infinite 
time)  linear  quadratic  regulator  is  natural.  The  main  advantage  for  considering 
the  infinite  time  case  is,  of  course,  that  the  optimal  control  gain  is  time 
inarlant.  Implementation  of  the  optimal  feedback  control  will  generally  require 
an  estimate  of  the  system  state.  For  the  stochastic  problem,  the  optimal  state 
estimate  is  provided  by  an  infinite  dimensional  Kalman  filter.  However,  it  is 
at  this  point  where  a  singlficant  gap  exists  between  theory  and  practice  in  that 
implementation  of  a  Infinite  dimensional  filter  is  not  generally  possible. 

The  most  popular  solution  to  this  dilemma  at  the  present  time  seems  to  be 
reduced  order  modeling  of  the  system,  see  for  example  [1],[2],[3],  and  [4]. 

Usually,  a  modal  representation  of  the  system  is  assumed.  The  basic  idea  is  then 
to  evaluate  various  computable  criteria  representing  the  significance  of  various 
modes  and  then  select  some  finite  subset  of  modes  to  represent  the  dominant 
dynamics  of  the  system.  A  linear  quadratic  regulator  (controller/estimator  com¬ 
bination)  for  this  finite  dimensional  model  is  then  designed  using  standard 
methods.  The  stability  of  the  system  using  the  resulting  compensator  is  usually 
checked  by  computing  the  eigenvalues  of  the  closed  loop  system  using  a  high  order 
19 truth  model”  to  represent  the  LSS  dynamics,  where  the  truth  model  is  of  finite 
order,  but  of  much  higher  order  than  the  compensator. 

1  Research  supported  in  part  under  Grant  No.  78-3550,  AFOSR,  USAF,  Applied  Math 
Division. 
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There  are  two  basic  problems  with  the  above  approach.  The  first  is  that  the 
approximation  is  performed  upon  the  system  rather  than  the  compensator.  In 
effect,  the  truncation  of  an  Infinite  number  of  modes  is  tantamount  to  "approxi¬ 
mating"  an  unbounded,  infinite  dimensional  operator  by  one  that  is  bounded  and 
finite  dimensional.  Thus,  the  infinite  dimensionality  of  the  problem  is  never 
used  in  designing  the  compensator.  The  second  problem  is  that  the  stability  of 
the  closed  loop  system  is  also  checked  against  a  finite  order  model  and  the 
closed  loop  stability  of  the  actual  infinite  dimensional  system  is  inferred 
essentially  on  faith. 

To  avoid  these  pitfalls  essentially  requires  that  we  have  a  partial  differ¬ 
ential  equation  (PDE)  to  represent  the  system.  For  complicated  structures  this 
may  be  asking  too  much.  Hence,  the  method  of  reduced  order  modeling  may  be  as 
viable  an  approach  as  any  in  those  cases.  However,  we  can,  and  should,  study  the 
problem  of  approximating  the  optimal  compensator  for  appropriate  systems  where 
the  PDE  is  known,  since  we  may  gain  insight  into  the  approximation  problem  that 
would  apply  to  more  general  and  more  complicated  systems. 

This  paper  considers  the  compensator  design  and  approximation  for  a  class  of 
systems  representative  of  large  space  structures.  An  analytical  solution  for  the 
compensator  is  obtained  using  the  infinite  dimensional  model  of  the  system.  A 
finite  dimensional  approx imat loti  of  the  compensator  is  then  derived  based  upon 
analysis  of  the  infinite  dimensional  compensator. 


2.0  SYSTEM  MODEL 

Ve  will  consider  the  following  system 

Mw(t)  +  A  w(t)  -  bu(t)  +  bn.(t)  (2.1) 

o  a 

y(t)  -  [b,w(t)]  +  e(t)  (2.2) 

(2.1)  is  the  inhomogeneous  equation  of  motion  of  an  undamped  oscillator  and 

(2.2)  is  the  system  measurement  equation.  v(t)  is  an  element  of  a  separable 
Hilbert  space  H  and  represents  the  small  displacements  of  the  system  (transla¬ 
tion  and  rotational)  relative  to  its  equilibrium  position.  The  operator  M  con¬ 
tains  the  mass  and  Inertia  properties  of  the  system.  M  maps  H  to  H,  is  linear 
bounded,  self-adjoint,  and  positive  definite.  AQ  represents  the  stiffness  of 

,  the  structure  and  maps  P(AQ),  the  domain  of  AQ  which  is  dense  in  H,  to  H.  A0  is 
linear,  self-adjoint,  closed  (or  can  be  closed)  and  generally  unbounded.  Ve 
assume  there  exists  e>0  such  that 

IAQ  w,w]  _>  E  IMI2  ,  W  E  t>( AQ) 

The  spatial  domain  of  AQ,  ft.  Is  bounded  and,  thus,  the  resolvent  of  A0  Is  com¬ 
pact  for  each  X  in  the  resolvent  set  of  A0.  Since  A0  is  closed  and  has  a  com¬ 
pact  resolvent,  its  eigenvalues  are  isolated  (countable),  have  finite  multipli¬ 
cities,  and  have  infinity  as  the  only  limit  point  ([5],  p.187).  In  addition, 
the  modes  k-1,2...}  of  AQ  are  orthogonal  and  form  a  basis  in  H. 

A  single  white  noise  disturbance  n^(t)  is  present  with  a  spatial  force  dis¬ 
tribution  defined  by  b  e  H,  The  support  of  b  is  assumed  to  be  small  compared  to 
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the  measure  of  ft,  so  that  bn<j(t)  represents  a  physical  realization  of  a  point 
disturbance.  A  single  control  u(t)  and  rate  sensor  y(t)  are  colocated  at  the 
sight  of  the  disturbance.  The  measurement  error  e(t)  is  a  white  noise  process 
uncorrelated  with  n^(t)  .  u(t),  n<j(t),  y(t),  and  e(t)  are  all  elements  of  R  for 
each  t  e  [O,«0.  Note  also  that  the  control,  measurement,  and  disturbance  are 
all  compact,  since  they  are  each  finite  in  number. 

In  order  to  utilize  semigroup  theory,  (2.1)  and  (2.2)  will  be  put  in  first 
order  form.  Let  xj(t)  ■  w(t)  and  x2(t)  «  w(t)  then,  the  system  (2.1)  and  (2.2) 
can  be  written  as 


where 


x(t)  *  Ax(t)  +  Bu(t)  +  Fnd(t) 
y(t)  -  Cx(t)  +  e(t) 


(2.3) 
(2. A) 


Note  that  P(A)  *  V[A0)  x  V,  where  V  *  {w  |A^  w  e  H).  Note  also  that  A  will  have 
a  compact  resolvent  since  AQ  does. 


As  a  function  space,  it  is  natural  to  use  the  energy  space  E  -  V  x  H  defined 
by  the  inner  product 

[xt*]g  *  [A^  Xj,  A^  z x]  +  [Mx2,  z2]  ,  x,z  e  E  (2.5) 


It  is  easy  to  verify  that  A  -  -A*  under  the  inner  product  (2.5).  Hence, 
ty  [6,  Corollary  A. 3.1],  A  generates  a  norm  preserving,  strongly  continuous 
group,  T(t) 

The  semigroup  can  be  represented  using  the  inodes  of  M”*A0.  The  modes 
are  orthogonal  under  the  inner  product  [M*,*]  and  are  complete  in  H.  Using  the 
tit's,  a  set  of  basis  vectors  in  E  can  be  defined  by  *  l^k»0]^t  +2k  “  [0*$klT» 
k-1,2,...  The  set  { ^2k*  is  complete  in  E  and  orthogonal  under 

\  !•••]£•  Thus,  any  x  e  E  can  be  written  as 


*  '  j£,*lk*lk  +  x2k*2k 

where  (4lk>x]£  -  xlk,  [+2k*X^E  “  x2k‘ 

A  representation  for  the  semigroup,  can  be  obtained  from  the  homogeneous 
solution  to  (2.3).  Writing  (2.3)  as 

*  -  “  £  <ilk‘X2k)  *lk  +  (x2k  +  \  *1^  *2k  “  0 


3 


we  see  that  Xj^  and  solve 


[:j  ■  U  :i  d 


for  each  k.  Thus  the  semigroup  can  be  expressed  as 
x(t)  -  T(t)  x(0) 


^  U  X2k(0) 

kTi  |Xlk(0)  C08“k *  +^“ 


I  ■  ’* .  +  rC  .tier)/ 
;  .  '  *  L  i  \b  • 1  •  Vv 


1 -  8lnu),  1 1 

“k  k  I  lk 


(2.6) 


+  I-wk  xlk(°)  8inwkt  +  *2k*°)  COSwkt5  ^2k[ 


3.0  STOCHASTIC  CONTROL  PROBLEM 

Given  the  system  (2.3)  and  (2. A),  we  want  to  find  a  control  u(t)  to  minimize 

1  (  -T  A  A  A  I 


J[u]  -  lim  7  E  |f  Iy2  |B*x(t)J2  +  Ju(t)|2]  dt 
T-h»  to  » 


(3.1) 


where  y>0  and  E{*}  denotes  expected  value.  Let  the  noise  covariances  be  given 
by  E  n^(t)  nj(t)  -  o|  I  and  E  e(t)  e*(t)  -  c*.  For  (A,B)  controllable*,  the 
minimum  of  (3.1)  is  attained  for 


u  (t)  -  -B*  P„  x(t) 

o  c 


where  x(t)  satisfies 

e 

x(t)  -  Ax(t)  +  Buo(t)  +  P^*  [y(t)  +  Bx(t)],  x(0)  -  0 

and  ?c  and  P^  satisfy,  respectively. 


lPcx,Ax]  +  [Ax,Pcx]  +  y  [B*x,  B*x] 


-  [B*  P  x,  B*  P  x]  -  0,  x  e  V(A ) 
c  c 


(3.2) 


(3.3) 


(3.4) 


[PjX,A*xJ  +  [A*x,Pfx]  +  Oj  fB*x,B*x] 

-  [B*  Pfx,  1*  Pfx)  -  0,  x  e  V(A*) 


*  Approximately  controllable,  see  [6,  Theorem  4.9.2). 
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(3.5) 


It  is  easily  verified  by  direct  substitution  that  under  the  energy  inner  product 

fc  -  yl  .nd  ff  -£l 

o 

thus  (3.2)  and  (3.3)  become 

uQ(t)  -  -  yB*  x(t)  (3.6) 

•  o 

x(t)  -  Ax(t)  +  BuQ(t)  +  — B[y(t)  -B*x(t)],  x(0)  -0  (3.7) 

d 

The  corresponding  minimal  cost  is 

Remark  3.1.  The  above  solution  requires  that  the  steady-state  Riccati  equa¬ 
tions  (3.4)  and  (3.5)  have  unique  positive  definite  solutions  and  that  x(t) 
be  asymptotically  stationary.  Regarding  (3.4)  and  (3.5),  since  we  have  ana¬ 
lytical  solutions,  the  only  concern  is  uniqueness  and  this  is  guarenteed  by 
(A,B)  controllability  [8].  The  asymptotic  stationarity  of  x(t)  follows  easily 
through  application  of  (6,  Theorem  6.7.1]. 

Remark  3.2.  For  morm  general  problems,  where  for  example  we  replace  B  in  (3.1) 
with  another  operator  R,  which  is  Hilbert- Schmidt,  exponential  stablllzabllity 
of  the  system  is  sufficient  to  guarantee  solutions  to  the  infinite  time,  sto- 
chastic  regulator  (ITSR)  problem  [6,  Theorem  6.9.1].  However,  when  compact 
controls  are  employed,  as  is  the  case  here,  the  system  cannot  be  given  a  uni- 
form  exponential  decay  rc.te*  unless  the  open  loop  system  is  already  uniformly 
exponentially  stable.  Note  that,  here,  the  open  loop  system  (2.3)  is  unitary, 
loS.,  l|T(t)||  “  1.  Hence,  the  closed  loop  system  can  only  be  strongly  stable** ***. 

Bamark  3.2.  The  existence  of  solutions  to  the  ITSR  problem  for  systems  that  are 
only  strongly  stabilizable  is  largely  an  open  problem.  Recent  results  concern¬ 
ing  strong  stability  and  the  steady— state  Riccati  equation  are  given  in  [8],  but 
the  sufficient  conditions  for  existence  are  quite  strong,  requiring  in  particu¬ 
lar  ||t*x||  <_  Mi  ||B*x ||  for  the  control  problem  and  ||F*x||  <_  M2||Cx||  for  the  filter- 
'  problem  mere  Mj  and  M2  are  constants?**  One  possible  physical  interpretation 
of  these  requirensnts  is  that  we  only  attach  a  performance  penalty  to  those 
points  where  we  have  located  actuators  and  that  we  locate  a  sensor  at  each  dis¬ 
turbance  source. 


*  A  semigroup  T(t)  is  uniformly  exponentially  stable  if  3  «>0,  M>1  3 
|Jr(t)  ||  <.Me“»t,  t >0. 

**  A  semigroup  is  strongly  stable  if  ||T(t)x||  -►  0  as  t-*“,  xeE. 

***R  is  the  state  weighting  operator,  F  is  the  disturbance  input  operator,  and  C 
is  the  observation  operator. 
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4 . 0  OPTIMAL  COMPENSATOR 


Let  T(t)  be  the  CQ  semigroup  generated  by  A.  Then,  the  solution  to  (3.7)  is 
given  by 


ct  a  I 

c(t)  -  J  T(t-x)  3u(t)  +  ^B[y(x)  -  B*  x(x)]|  dx 
o  d  ' 


(4.1) 


Using  (3.6),  we  can  rewrite  (4.1)  as 

t 


[yy (T)  -  b*  x(x)](  dx 


u(t)  ■=  -  f  B*  T(t-t)  B  jyu(x)  +  ^ 

*  o  ’  d 

-J  B*  T(t-x)  B  +  y(T>|  dx 


Therefore 


t  ^ 

u(t)  +  k  f  B*  T(t-x)  Bu(t)  dT  -  -  k  f  B*  T(t-x)  By(x)  dx  (4.2) 

u  y 


where  ku  -  y  +  oa/o&,  k  *  yoB/od.  Observe  that  (4.2)  is  a  Volterra  equation. 
Define  the  Volterra  operator  L  as 


t 

if  *  h  ;  |  B*  T(t-r)  Bf(r)  dx  m  h(t) 

o 

Then  (4.2)  can  be  written  as 


u  +  k  Lu  “  -  k  Ly 

u  y 


(4.3) 


Solving  for  u  we  obtain  the  optimal  compensator  in  abstract  form 


uc  -  -  ky  [1  +  ku  L]~l  Ly 


-  -  ky  J1  +  K]  Ly 


(4.4) 


'where  K  is  also  a  Volterra  operator  of  the  form 

Kf  -  h  ;  J  K(t,x)  f(x)  dt  -  h 
o 

Prom  [6,  p.102-103],  K(t,x)  can  be  computed  iteratively  and  is  given  by 

K(t,x)  -  E  *„(*’*> 

where 


n-1 


(t,x)  -  gn(t-x)  -  J  g(t-d)  gn_j (o-t)  da 


(4.5) 


(4.6) 
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gjtt-T) 


g(t-x) 


B*  T(t-x)  B 


(4.7) 


Noting  that 

Ln f  *  j  gn(t-T)  f(x)  dx 
o 

(4.4)  can  also  be  written  as 

u  *  -  k  -  k  +  k^  L2  -  k^  +  . . . 1  L  y 

c  y  (  u  u  u  ) 

-  -  ky  jr  (-fcu)n_1  tn|y  (4-8) 

Since  B*  T(t-x)  B  is  uniformly  continuous  and  (4.3)  is  a  Volterra  equation,  the 
iteration  defined  by  (4.5)-(4.7)  and,  hence,  the  series  (4.8)  converge  for 
te[0,Tl *T<t,T<  •,  [6]*  The  series  (4.8)  is  sometimes  called  a  Neumann  series, 

[9].  Thus,  the  optimal  compensator  is  of  the  form 

u(t)  -  -  k  J  go(t-x)  y(x)  dx  (4.9) 

y  *  o 

where  the  kernel  gQ(t-x)  is  given  by  the  bracketed  term  in  (4.8).  Note  that 
(4.8)  utilizes  the  infinite  dimension  system  model  via  (4.7)  which  is  given  in 
terms  of  the  open  loop  semigroup  operator,  T(t),  and  we  have  a  representation 
for  T(t),  (2.6),  using  the  modes  of  M~^A0.  In  particular,  expanding  (4.7)  on  the 
orthonormal  basis  for  E,  as  defined  ln  Section  2,  gives 

g(t-x)  ■  B*  T(t-x)  B  -  53  bv  COSU.  (t-x) 

k-1 


where  bfc  -  [M(W1b),  $2kl  ”  lb*  ♦fcl* 

Properties  of  g(o) 

i)  g(o)  •*>  0  as  a  -*■  • 
il)  g(o)  is  uniformly  continuous 
ill)  g(o)  is  compact  for  each  t  _>  0. 

Condition  1)  follows  since  g(o)  “  B*  S(o)  B  and  S(o)  is  the  CQ  semigroup 
generated  by  A  -  k^BB*  which  is  strongly  stable.  Note  that  S(o)  is  also  a  con¬ 
traction,  l.c.,  ||S(0)||<  1.  The  second  condition  results  from  the  strong  con¬ 
tinuity  of  8(ff).  For  condition  ill),  g (a)  is  compact  since  B  is  compact  and 
S(o)  is  a  bounded  operator. 
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5.0  SMALL  GAIN  APPROXIMATION  TO  OPTIMAL  COMPENSATOR 

Using  the  procedure  outlined  in  Section  4.0,  four  complete  interations  and 
a  partial  fifth  iteration  were  performed.  This  yielded  an  approximation  for  the 
kernal  of  the  optimal  filter  given  by 


gn(cr)  *=  cosed,  o 

k-1 


[i  -  cck«  +  |(cck»)2  -  jr<ccto)3  +  1 


4  sck 


] 


-  £  {ck  bk  [>  -  W- +  Kk»>  -  3T<W»:  ] 


00  CD 


-  E  E 

k-1  j-1 

k*J 


bv  Sinai,  a  r 

kj  l1  " 


w  +  Kk")2i  +  Qi>c  . 


■] 


sinm^a 


[’-V'TV’2]  +  k«  * 


a) 


5.1) 


where 


i  I  "  R.  D, 

ck  u  k 


k„  b^/2  .  *  u>k  bj/(“k  -  Uj)  ,  a  -  t  -  r 


2  1 

The  remainder  term  0^(1^  ;  55—  ;  a)  is  given  in  the  appendix  and  is  composed 
of  terms  of  second  order  or  greater  In  k^. 

It  Is  apparent  that  the  terms  in  the  square  brackets  in  (5.1)  are  the  first 
few  terms  of  the  Taylor  series  expansion  for  exp  (-£ckt) .  Thus,  the  limiting 
form  of  the  kernel  is 


where 


*o(o)  “  *o(a)  +  ku  0(ku;  7"  ;  o) 

k 


*o(o)  "  S.  8lk(a)  *  ku  I*  E  «2kj(o> 


k-1 


(5.2) 


(5.3) 


k*j 


2 


*2kJ<,’>  *  °kj  k.  • 


2  -{.k« 


[einut.  o  "] 

~~k»-'ck  — J 

(=?=)•*»>; 


(5.5) 


Kot*  that  tha  coaponant  flltars  gjfc(a)  have  second  order  dynamics  while  the  fil- 
ters  S2kj(°)  have  fourth  order  dynamics.  Also,  all  the  parameters  in  gik(o)  are 
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related  to  the  kth  open  loop  mode,  while  *s  re^atec^  to  t^ie  cross-coupling 

between  the  kth  and  jth  modes.  Due  to  the  presence  of  the  factors  l/a)^  and  1/wj 
in  g2kj(°)»  the  cross-coupling  terms  are  attenuated  relative  to  the  "diagonal” 
terms  glk(t~T)  as  wk  and  «*>j  become  large.  Thus,  at  high  frequencies  the  dominant 
filter  dynamics  are  due  to  the  diagonal  terms  g^(a).  Since  the  iteration  con¬ 
verges  and  since  gQ(a)  and  gQ(a)  are  continuous  in  a,  k£  0(ku;  1  /oj^jq)  is 
continuous  in  a.  Also  k^  1 0(ku;  l/w^ja)  |  -*  0  as  ku  •>  0  for  each  ae[0,T].  Therefore, 
for  any  e  >0,  we  can  find  sufficiently  small  such  that 

max  |k^  0(k  ;  —  ;o)  I  <  e 
0<t<T  k 


The  argument  is  used  to  indicate  that  [OC^; l/t%;a)  |  0  as  i%  ^  Therefore, 

assuming  ky  is  small,  the  approximate  optimal  control  is  given  by 


-k 


£ 

k=>l 


V- 

I 


glk(t-T) 


y  (t)  dT  +  k  k  £ 


C  £  J 

=  1  j*=l  Jc 


82kj(t"T^y(T)dT 


(5.6) 


k?*j 


Since  ky  and  1^  are  both  small,  as  crude  first  approximation,  the  double  sum¬ 
mation  in  (5.6)  could  be  neglected.  In  this  case,  the  filter  takes  a  particularly 
simple  form,  namely, 


u(t)  =  ^(t) 

k=l 


(5.7) 


where 


Uk(t) 


8lktt-T)  y(T) 


dr 


Here,  the  basic  structure  of  the  compensator  is  an  infinite  bank  of  second  order 
filters,  each  operating  on  the  sensor  output  and  the  control  is  simply  the  sum 
of  their  outputs*  The  component  filters  gik(t-T)  in  this  case  have  a  rather 
special  property.  Letting  Gik(s)  denote  the  Laplace  Transform  of  gjjc(t-T),  we 
have 


Gik(,) 


<• 


■2  +  2£ck* +  4 +  \ 


Such  a  filter  is  strictly  positive  real,  [111,  [12].  Note  that  the  phase  shift 
of  the  filter  never  exceeds  ±ir/2.  According  to  [13),  strictly  positive  real  fil¬ 
ters  have  robust  stability  properties  in  the  sense  that  they  will  stabilize  any 
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positive  real  system  (in  an  input/output  sense).  The  system  considered 
here  is  positive  real  since  it  is  an  undamped  oscillator.  Note  that  any  finite 
sum  of  positive  real  filters  is  also  positive  real.  Therefore,  we  can  truncate 
the  summation  (5.7)  at  any  finite  number  of  terms  and  the  compensator  will  still 
stabilize  the  system.  It  would  be  expected  that  the  more  terms  retained  in  the 
summation  the  closer  the  resulting  performance  would  be  to  the  optimal  performatce, 


Returning  to  the  filter  defined  by  (5.6)  another  approximation  will  be 
obtained.  First,  consider  the  Laplace  Transform  of  ku  g2kj(l)  given  by 


k  G0..(s) 
u  2k  j 


r  2 

hi  njk 

(s  +  eck)2  +  w2 

(s  +  Ccj)2  +  o.2_ 

Using  the  definitions  of  £ck*  % j »  etc,  this  becomes 


ku  G2kj(s) 


bk  b1  {(h>k  ~  “ft  *  +  2ku(uk  b1  ~  hb  +  ku((Jk  b1  "  “l  bl)/4} 


<V 


Dropping  all  terms  second  order  or  greater  in  k^  gives 


kuG2kj(s)  “  ku  ,  2 


.2  .2  2 

bk  b1  S 


2TS  kuClk<8)  Glj(8)  (5‘8) 


(®‘+  kub^8+  u£)  (s2+  kub2s+  up 

2  2 

where  Glk(s)  is  the  Laplace  Transform  of  glfc  (t)  (neglecting  £ck,  ?cs  in  the  denom¬ 
inator)  .  -Thus  the  Laplace  Transform  of  the  compensator  is 


Go<8>  "  "kv  £ 


b£8 


.2  ,2  2 
b  b  s 
k  P1 


»  «  “  S  j?i  <.2*  b 2  +  J)  (.2*  b2,  ♦ A 

uk  k  kltj  uk  k  Uj  j 


(5.9) 


In  the  time  domain,  the  approximation  to  the  optimal  control  is  now  given  by 

G(t)  -  -ky  £)  *lk(t_T)  y(x)  dT  +  kyku  £  J  jj  glk<t_o)  8ij(o_T>d®jy(T)dT 


krj  (5.10) 

This  representation  is  Interesting  in  that  the  basic  building  blocks  of  the  fil¬ 
ter  are  the  positive  real  component  filter  gik(t-t).  Note,  however,  that  the  con¬ 
volutions 
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83kj(t,T)  "  J  8lk<t_0>  ®ij (°~T)  da 
or,  equivalently,  the  products 


emw  '  °ik<8)  Gu<8) 

are  not  positive  real*  However,  since  the  diagonal  terms  dominate  at  high  fre¬ 
quencies,  the  high  frequency  component  filters  approach  a  positive  real  form. 
Since  we  would  expect  a  close  (finite  dimensional)  approximation  to  the  optimal 
filter  to  stabilize  the  system,  it  is  interesting  that  the  limiting  structure 
is  positive  real,  i.e.,  a  structure  that  is  compatible  with  truncation. 

6.0  CONCLUSIONS 

It  has  been  shown  that  the  optimal  compensator  can  be  represented  in  terms 
of  a  kernal  of  an  Integral  operator.  The  results  demonstrate  that  the  integral 
representation  provides  considerable  insight  into  the  structure  of  the  compen¬ 
sator.  This  Is  due  to  the  compensator  being  compact  and,  hence,  something  that 
is  Inherently  approxlmatible  by  a  finite  dimensional  operator,  i.e.,  the  compen¬ 
sator  parameters  go  to  zero  as  k  -*■  •.  For  the  system  considered  here,  the  opti¬ 
mal  filter  can  be  constructed  of  basic  building  blocks  which  are  second  order 
positive  real  filters.  The  first  (low  gain)  approximation  is  a  diagonal  array 
of  these  component  filters  and  there  is  a  one  to  one  correspondence  with  the 
open  loop  modes.  The  second  (low  gain)  approximation  adds  in  the  cross-coupling 
effects  in  the  form  of  fourth  order  filters  which  are  convolutions  between  the 
various  component  filters  taken  two  at  a  time.  At  high  frequencies,  the  dynamics 
of  the  conq>ensator  approach  the  diagonal  form  of  the  first  approximation  which 
should  have  advantages  when  the  filter  dynamics  are  truncated  In  order  to  obtain 
a  finite  dimensional  approximation  for  the  compensator.  The  results  also  seem 
to  suggest*  that  higher  gains  will  require  Increasingly  higher  order  convolutions 
of  the  component  filters. 
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APPENDIX 

The  following  defines  O^kyit),  the  remainder  term  for  the  optimal  compen¬ 
sator,  after  four  complete  iterations  and  part  of  a  fifth  Iteration: 
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w  -  E 


b^t  sinw^t 


k-1  k 


«k1  bk8lnwkt  /  tbk  \  nikb1sitlw1t 

♦  EE  ('i(k-3>+  -r(1  -*<*'>)  *  1  ■ 1 


ut  + 

*ck 


Kk‘)2] 


(cjU.kH^Ci-s^t)) 


♦£F  ❖*  l(c2(k,j)+kuc3(k,j)t)  coso^t 
krj 


+  (c,(j,k)Anc,(jlk)t)  cosu.t]  -ELE  bSJbJ  I(c.(k,j,i)4k  c s(k,J,i)t) 

2  u  3  3  k*J*i  k  3  1  4  u  3  . 

.  cosw^t  +  (c4(j,k,i)+kuc5(j,k,l)t)  cos<i)jt  +  (c^i,  j  ,k)+kuc5(i,j  ,k))  cosujtlj 

♦  k A  i  .  ece  (Vk.,.» 

u  (k-l  ®“k  12a>2  k^j^i  k  3  1  6  ^ 

+  c,(j,k,i)  siau  t4c,(l,j,k)sinw.tl  - EEEEb?b^b?b 2  tc7(k,j,i,l)  sinu^t 
6  j  6  i  k*j*i*l  k  3  1  1  1 


+  C7(j,k,l,£)  sinu^t  +  c7(l,j,k,£)  sinWjt  +  c7(£,j,i»k)  slnu^t] 


where 


2  \  4  2  “X  3bX“jl  2  K-X 

c.(k,j)  — rV?  bX  +  2  —  HH-  -~r~? 

<V“J>  [3k  «2  0,2  J  sh  L  -k 

c2(k,J)  - - 2~3  j  -2-  <w?b2  -  «£b2) 

2  2(u£-«2)2  3k  k3 


‘J  *•£ 


c3(k»j)  "  77 
(tV 


1 _  \.(  2  2  4  “khfchjl  3bk 

-»>ll  v  i  ^  2  J'^2 


c4(k.3)  -  2  r .  2  2.  •  cs<k»J»i3  ""TTJTT 
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c6(k,j,i) 


Cj(k|  j  »i»l) 


,  2  2W  2  2. 


,  3,2  2{“k(“i'Hil?)‘  2Vj“i)b£  2 

Hbi + - rH; — J - +  “kbk 


<V# 


3 

“k 


.  2  2W  2  2W  2  2. 

<ui"“k)  (tl)j‘“k)  (ut’°)k) 


*ck  *  kubk/2  *  Rkj  “ 
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